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Abstract. For a fairly general class of two-dimensional tiling substitutions, 
we prove that if the length expansion (3 is a Pisot number, then the tilings de- 
fined by the substitution must be locally finite. We also give a simple example 
of a two-dimensional substitution on rectangular tiles, with a non-Pisot length 
expansion /3, such that no tiling admitted by the substitution is locally finite. 
The proofs of both results are effectively one-dimensional and involve the idea 
of a certain type of generalized /3-transformation. 



1. Introduction 

This paper was motivated by the study of simple tiling substitutions such as 
the substitution S shown in Figure ^ To the left of the arrows, four prototiles 



T 



I 1 



I'D- 



Figure 1. An inflate-and-subdivide rule S. 



are shown along with their dimensions. The larger dimension, (3 = (1 + \/13)/2 « 
2.3027756377, is a root of the polynomial p(x) — x 2 — x — 3. The patches to the 
right of the arrows show the application of the substitution S to each prototilc. 
The substitution S can be applied repeatedly, to produce arbitrarily large patches 
of tiles. The first few applications, starting with the large square prototile, are 
shown in Figure Our interest in S arises because its length expansion (3 is not 
a Pisot number. An algebraic integer (3 is called a Pisot number if j3 > 1, and all 
its Galois conjugates (3 1 satisfy |/3'| < 1. Here, (3 is not Pisot since its conjugate is 
/3' = (1 - Vl3)/2 « -1.3027756377. 

It is known that whether or not the linear expansion of a tiling substitution S 
is a Pisot number has implications for the tiling dynamical system that it defines. 
(We make the distinction between the linear expansion and the length expansion 
of a substitution tiling precise in Section 0J) Solomyak proved in [5] that a one- 
or two-dimensional self-similar tiling dynamical system with real linear expansion 
(3 is weakly mixing if and only if f3 is not Pisot. A similar result holds for a 
two-dimensional self-similar tiling dynamical system with complex linear expansion 
A G C\M. Solomyak showed that such a dynamical system is weakly mixing if and 
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only if A is not a "complex Pisot" number. However, in the two-dimensional case the 
proofs of these results, as well as many similar results, use the additional assumption 
of "local finiteness" (a property that is automatic for the one-dimensional case). A 
tiling is defined to be locally finite 1 if, up to rigid motions, it contains only finitely 
many two-tile configurations. 

For substitution tilings, local finiteness follows from "two-tile closure property" (i . 
On the other hand, Danzer has conjectured that non-locally finite tilings may 
in fact be generic for certain classes of substitutions. We will show that tilings 
admitted by the substitution in Figure [2 fail to be locally finite. It remains unclear 
the extent to which results like those in [9] can be extended to the non-locally finite 
case. 

Although some examples of non-locally finite tilings have already appeared (see 
PQ, see also pj), our example has the benefit of simplicity. It belongs to a family 
[2] of tilings that can be created quite easily from any one-dimensional substitu- 
tion. Moreover, our proof is based on a one-dimensional method that may be of 
independent interest. 

To see how the lack of local finiteness appears from Figure consider Figure 
which shows three iterations of the large square tile. In the first iteration, we have 
circled a point where there is, in Danzer's terminology pp, a "misfit situation". 
Under the inflate-and-subdivide rule S the tiles are inflated by the linear map 



TX 



TX 



TX 



TX 



TX 



TX 



LT 



LT 



LT 



Figure 2. Three iterations of the larger square tile. 

(31; we can follow the circled point's location under this linear map. Following 
the circled point through two more iterations, we see that the misfit it creates 
moves around relative to the tile below it. We will prove that this point occupies 
infinitely many locations with respect to the tile below it. Therefore the infinite 
tiling generated by this process is not locally finite. 

Now one can begin to see why the arguments are one-dimensional. Inspection of 
Figure^shows that there are two kinds of upper edges: a long one and a short one, 
which we denote I and s respectively. The reader can verify that a long upper edge 

^This property is also sometimes also called "finite local complexity". In the theory Delone 
sets, however, "local finiteness" usually refers to a weaker property. 
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is substituted by a long followed by three short edges, (I — > Isss), and a short upper 
edge is substituted by a long one, (s — > I). The circled point in Figure[21is always on 
the corner of two of the tiles above it, but it is (we will show) always in the interior of 
an upper edge. Thus we restrict our attention to the one-dimensional substitution 
on the upper edges to prove our result. In fact, following the misfit through these 
iterates is analogous to applying a suitably defined "generalized /^-transformation" , 
which we discuss in the next section. 

The organization of this paper is as follows. We first focus on the one-dimensional 
theory, which is itself of independent interest. In Section |2 we give the one- 
dimensional definition of a substitution tiling and a generalized /3-transformation, 
and we show how they are related. In Section we give conditions that guarantee 
finite or infinite orbits of the generalized /3-transformations and relate this result 
to one-dimensional tilings. In Section 01 we apply this machinery to tilings of M 2 . 

2. Definitions 

2.1. One-dimensional tiling substitutions. Let A = {1, ...,m} and let A* 

denote the set of all (nonempty) finite words in A. A substitution is mapping 
o~ '. A — > A*. Given a substitution a, the structure matrix is the d x d matrix M 
with entries Mi j equal to the number of occurrences of the symbol j in the word 
o~(i). We require that M has a real eigenvalue j3 > 1 such that the right eigenvector 
h = (hi, hi, . . . , h m ) has strictly positive entries. It is worth noting that Perron- 
Frobenius theory ensures this in two special cases: irreducibility and primitivity. 
We say M is irreducible if for each i,j there is an n for which (M n )ij > 0. If a 
single n works for all i,j then we say M is primitive. By the Perron-Frobenius 
theorem, if M is irreducible then there exists an r > for which h G (rQ(f3)) m . 

Given a substitution a, we construct a one-dimensional tiling substitution S and 
a corresponding tiling r of the line. The prototiles V = {Pi, ■ ■ ■ , P m } are intervals 
Pj = [0, hj) with lengths given by the entries of h. A tile is a translation of a 
prototile T = Pj + s = [s, s + hj). A tiling patch ir is a finite sequence of tiles laid 
end to end: 

(1) = {P n + S l) P J2 + *2i • • • . Pj n + s n) 

where 

Si+i — Si + hji- 

Equivalently, in a tiling patch, the tiles are disjoint and the support, defined 

(2) supp(Tr) := U? =0 P i( + m, 

is an interval. A tiling r is an infinite set of disjoint tiles whose support is M. 

We call s the left endpoint of the tiling patch tt, and if supp(7r) = [s, t), we call 
t the right endpoint. We can specify a tiling patch tt by giving its left end point s 
and the word w = j\ ]2 ■ ■ ■ j n that specifies the prototiles that occur in its sequence. 
In particular, the patch in Q is given by tt — 7r(s,w). 

Now we can define the tiling substitution S corresponding to a. It is the mapping 
from tilings to tiling patches defined by 

S(P j + s)=ir((3s,o-(j)). 

Clearly we have supp(S'(T)) = /3supp(T), and S(T + s) = S{T) + (3s. This allows 
us to extend S to a mapping on tiling patches, and even to tilings. A tiling r is said 
to be ^-admissible if every patch tt C t is equivalent to a patch of tt' C S n P for 
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some P E V, n € N. The set of all S'-admissible tilings (which is always nonempty) 
is called a tiling space, and is written Xs- 

Example 1. Consider the substitution I — * Isss, s — ► I from the introduction. The 
structure matrix is 



which has characteristic polynomial p(x) = x 2 — x — 3, and the Perron eigenvalue 
/3 = (1 + V13)/2. Taking h = (/3, 1) as a positive right eigenvector, the prototiles 
are P\ = [0, (3) and Pi = [0, 1). These tiles have the same lengths as the edges of 
the tiles in Figure Q We will continue our discussion of this example below. 

2.2. Generalized /3-transformations. Fix an interval E = [0, B) and a real num- 
ber (3 > 1. 

Definition 2.1. We say a mapping Fp : E — > E is a generalized /^-transformation 
if (i) Fp is piecewise linear, (ii) Fp is continuous from, the right and (Hi) F'^(x) = (3 
for almost every x E I. 

More explicitly, Fp is a generalized /3-transformation if there is a partition of E 
into I subintervals Ij = [ ), with endpoints = x a < X\ < ... < x\ — B, and 

real numbers yi, j/2, ■■■Vi with f3Ij — yj C E for j — 1, . . . , I, so that 

(3) Fp{ x ) — 0% — yj whenever x G Ij. 

One example of this is the classical /^-transformation Tp : [0, 1) — * [0, 1), defined 
Tp{x) = (3x mod 1 (see Renyi 0). A type of generalized /3-transformation, more 
general than what we consider here, was studied by Wilkinson JT]. He showed 
under several additional conditions that Fp has an absolutely continuous invariant 
measure and is weakly Bernoulli with respect to the partition into the intervals Ij . 

Definition 2.2. We call a generalized (3-transformation Fp algebraic if (3 is an 
algebraic number and there is an r € R for which 

x 0) ■ ■ ■ ,xi,y 1; . . . ,yi £ rQ(J3). 

2.3. Substitution tilings and generalized /3-transformations. Now let us 

think of the circled point in Figure as being the origin situated in the interior of 
the one-dimensional tile that is the upper edge of a two-dimensional tile (Figure OJ. 

E • ) 

E )E ) 

I 

f )E * * •Ht * ) 

Figure 3. The map S finds the central tile after the substitution is applied. 

Suppose that T = [s, t) is a tile with s < and t > 0, so that S T. We call 
such a tile a central tile, illustrated in Figure with T = [1 — j3, 1). It follows 
that G supp(S'(T)) = [/3s, (3t), and thus there is a unique central tile in S(T) = 
(Pj 1 + si,Pj 2 + s 2 , . . . , Pj n + s n ), i.e., € Pji + Si for some i. We can therefore 
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define a mapping S from the set of all central tiles to itself by taking S(T) to be 
the central tile in S(T). We are interested in the effects of iterating S on central 
tiles, but it is convenient to think of this process a little differently. 

First we consider the patch it := tt(0, 123 . . . m) = (E\, E2, . ■ ■ , E m ) where Ej = 
Jj + gj with go = and gj+i = gj + hj. We write the support of tt as E — [0, B). 
The mapping F : E — > E is defined as follows. For x G E, find j G A so that 
x G -Ej. Then Ej — x is a central tile, and S^-Ej — x) = E& — y for some unique 
k G A and t/ G M. We put y = F(x), noting that F(x) G Efc. 

We may write F explicitly as a generalized /3-transformation by partitioning 
each Ej into subintervals that will be sent by F onto whole tiles (the canonical 
subdivision suggested by S). We obtain a partition = xq < x\ < ... < xi = B 
into intervals 7i, such that for each If. there is a j such that F(Ik) = Ej. The 
restriction of F to Ik is the linear map taking Ik onto Ej . So for x G Ik we have 

(4) F(i)= ) 9(x-a!fc_i)+s J -_i; 

hence y k = -0Xu-i +9j-i- 

Example 2. The generalized /3-transformation for Example is defined on the 
interval E = [0, 1 + /3) so that £4 = [0,/3) and E 2 = [/?,/?+ 1). The partition of 
E given by the substitution is xq = 0, Xi = 1, £2 = 1 + l//3,#3 = 1 + 2//3, IE4 = 
1 + 3//3 = /3, and = 1 + (3. Since go — and gx = ft our translations are 
2/1 = 0, t/2 = 0, J/3 = 1, 2/4 = 2, and ys = 3 + 0. The transformation is as shown in 
Figure 0] 




1 B 1+B 

Figure 4. The map F : E — > E 1 . 

Since each interval in a tiling substitution expands by a factor of /3, we have the 
following lemma. 

Lemma 2.3. The mapping F given by a tiling substitution with expansion constant 
(3 is an algebraic generalized (3 -transformation whenever h G (r(Q>(/3)) m . 

This is certainly true whenever M is irreducible or primitive. 
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3. One-dimensional results 
Suppose (3 > 1 is a root of the irreducible monic polynomial p(x) = go + gix 



92% 



{l,/3,/3 2 ,.../3 d - 1 }oveT 
by A in Q d , where 



(5) 



G where d > 1. We can identify 



with 



via the basis 



Multiplication by /? in 



corresponds to multiplication 



.4 



1 - 



-go 

-Si 
-52 



is the companion matrix of p(x). 

To fully understand the action of A, we embed 



in C d (although if the eigen- 



values are all real, it suffices to embed in R d ). We define a surjective linear map 
<f> : C d -> C by 

<f>(w) = w + + w 2 /3 2 + ••• + u'tf-i/^ -1 . 
The restriction of <\> to Q d is the isomorphism between Q d and Q(f3). Moreover, for 
any w G C d , 

(6) = /30(u?). 

The Galois conjugates /3a, /?3, ■■■/3d of /? are eigenvalues of A with eigenvectors 
V2,V3, -.-Vd- By linearity ^(j4u,-) = 4>{(3jVj) — (3j4>{vj), but by 10 we also have 
<t>{Avj) — (3<t>(vj), and so = for all j = 2,3, ...,d. Let Vo be the subspace 

of C d spanned by V2, ■ ■ . , Vd, arid let V + be the (one-dimensional) eigenspace of (3. 
Since <fi is not the zero map, we can select the vector v G V + with <fi(v) — 1 as a 
basis vector. 

Consider fixed an algebraic generalized /3-transformation F : E — > E. Recalling 
the notation of Subsection 12 . 21 we have the partition of E given by = xq < x\ < 



< xi = B with Xj € 



; we note that by definition 0{xj 



)<B. 



denote the jth subinterval as Ij = [xj-i, Xj), and we have y\, j/2, ■■■Hi € 
that 13 1 j - yj C E. 

Since <f>(u) = for any u £ Vq, we know that </>(> 
Thus we can define 



We 
so 



u) € R if and only if c € 



W[ a}b ) = {w : <f>{w) G [a, b)} = {cv + u : u e V and c G [a, 6)} C C d 

to be the set of all points in C d that <fr maps onto the real interval [a, b). 

Let We denote the subset of C d corresponding to E and denote the subsets of 
C d that correspond to the intervals Ij as Wj — Wza the WjS form a partition of 
We- There are unique vectors w\, W2, Wd in Q d for which <j>{wi) = y,, and one 
can easily check that AWj — Wj C W^g. Thus we can define the map / : We — * Wb 



(7) 

and we see that 
(8) 



f(w) = Aw 



uij whenever w G Wy, 



(f> o f = F o cj). 



Because our interest lies with the mapping F. the theorem below is stated in 
terms of it, but the result is actually proved for / and then extended to F. The 
proof is mostly linear algebra and uses ideas similar to 0, [E] and ^H] ■ 
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Theorem 3.1. Let (3,(32, ■■■,(3d, E, and F be defined as above. 

(1) If \(3j\ < 1 for all j = 2, 3, d, then for any x £ Q(/3) n E the orbit of x 
under F is finite and therefore eventually periodic. 

(2) If \(3j\ > 1 for some j £ {2, 3, d}, then there exists an x £ Q(/3) H 13 such 
that the orbit of x under F is infinite. 

Proof. Since <j> o f = F o (j) and the map is a bijection when restricted to Q d it 
suffices to prove the results for /. 

Suppose | < 1 for all j = 2,3, ...,d. Let w £ Q d n We and uJi from Q be 
written in terms of the eigenvectors of A as 

d d 
w = cv + CiVi and iOj = + Wj^Vi. 

i=2 i=2 

For i = 2, . . . ,d choose 

n > fZTgj fora11 iei,2,...,d 

Define the trapping region IZtrap to be the set of all points in such that Cj < rj. 

We claim that IZtrap is invariant under the action of /. Suppose w £ TZ tr ap with 
such a choice of values, and find j £ 1, 2, I such that w £ W}. Then 

(9) /(iZf) = Aw- Wj = + y^jCjVi) - (wj,iv + y^Wj^Vj). 

i=2 i=2 

Since the component in the direction of v is unrelated to the definition of TZtrap 
we see that we must have |/3jCj — Wji\ < r< for i = 2,3, ...d. This is true since 
\ w j,i\ < (1 — an d |cj| < r,. This proves the invariance. 

Now note that there is some q £ Z for which uJ and the WjS are in |Z d . Since A is 

an integer matrix it leaves "^L d invariant and so / (j^> d H Wej C ^Z d r\WE- Choose 
the numbers for the trapping region so that to € TZtrap] clearly / ( -Z d D 1Z tra p ) !_ 



9 

iZ d R TZtrap- Since ^Z d n TZtrap is a finite set the orbit of w is finite and thus 
eventually periodic. 

To prove the second part of the theorem we must show there is an infinite escape 
region TZ esc C We for which orbits beginning in TZ esc tend to infinity. Without loss 
of generality we may assume that |/?2| > 1- Choose r 2 so that 

r2> uTT 1 T fora11 3 = ia,-,d. 

\P2\ - 1 

We define 

d 

TZ esc = {w = cv + CiVi C We such that \c<i\ > r 2 }. 

i=2 

To prove the result, we must show that for any w £ TZ esc H |Z d , the component 
of f(w) in the direction of v 2 exceeds that of w in modulus. Supposing again 
that w £ Wj, by the computation in equation JJjJ this amounts to showing that 
\(3 2 c 2 ~w^ 2 \ > |ca|. By choice of r 2 we have that \w 3t2 \ < (\(3 2 \~l)r 2 < (\{h\ ~ 1) j c 2 j 
and the result follows. 

□ 
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Corollary 3.2. Let S be a one- dimensional tiling substitution with expansion con- 
stant (3 > 1 and h G (rQ(/5)) m . Let T be a central tile, (i) Suppose (3 is a Pisot 
number. If the endpoints of T belong to Q(/3), then there exist K, L € N so that 
gNL+K (rp^ _ c;K(rp^ j or a n > 0. (ii) Suppose (3 is non-Pisot. Then there is a 
q € N and finite E C |Z[/3] smc/i that for any s G -Z[/3]\i? ; i/T = [— s, t) is central, 
then S K {T) ^ S L (T) for any K > L>1. 

4. Tiling Results 

In this section we study a class of self-similar substitution tilings in R 2 = C. Our 
first result is that if a substitution S has a Pisot length expansion, then the tilings 
it produces are locally finite. This solves a conjecture of Danzer pp for the type of 
tilings that we consider. In this section we also prove that tilings generated by the 
substitution from Figure ^ are not locally finite. 

We begin with some general definitions. A tile T is a closed topological disk 
DCC, together (possibly) with a label which is an element of some fixed finite set 
A. We say D = supp(T) is the support of T. Two tiles are called equivalent if they 
carry the same label and have congruent supports. 

We will fix a finite set V of inequivalent tiles, called prototiles, and assume all 
tiles T are equivalent to tiles from V . A packing is a collection tiles that intersect 
only on their boundaries, and the support of a packing is the union of the supports 
of its tiles 2 . A tiling patch it is a finite packing with a connected support. We 
denote the set of patches by V*. A tiling r is a packing with support C. 

Let A G C with |A| > 1. Suppose S : V -> V* is such that for each P G V, the 
supp(5(P)) is congruent 3 to A ■ supp(P) := A • P. We call S a self-similar tiling 
substitution (or simply a substitution when there is no danger of confusion). We 
call A the linear expansion of S, equating multiplication by A with a linear map on 
R 2 ; volumes must expand by a factor of AA. We call j3 = |A| the length expansion 
of S because it gives the expansion of one-dimensional objects. 

The substitution S extends in an obvious way to any tile T, and can thus be 
iterated to obtain infinitely large patches S n (P). Following the one-dimensional 
case, we say a tiling r is S-admissible if for every patch n G r is equivalent to a 
patch 7r' G S n (P) for some P G V and n > 0. We denote the set of S'-admissible 
tilings r by Xg. 

The following Lemma may be well-known, but we include a short proof for 
completeness. 

Lemma 4.1. Let V be a finite prototile set and let S be a tiling substitution on V . 
No prototile in V contains a nontrivial circular arc in its boundary. 

Proof. Suppose the boundary of a prototile P G V contains an arc of radius r > 
subtending an angle 9 > 0. The boundary of S n (P) contains arcs of radius (3 n r 
subtending the angle 9, and this implies that there are prototiles in V that have 
circular arcs of radius (3 n r in their boundaries. Since the prototile set is finite, we 
can find the maximum angle 9 n > subtended by the arc of radius P n r in any 

2 In the one dimensional case, it was convenient to take tiles as half-open and insist on dis- 
jointness. This also forced equivalence to preserve orientation, a restriction we do not impose 
here. 

■^In the event that there is more than one congruence between the supports S(P) and A • P, 
we will use the implicit parameterization of the boundaries to ensure uniqueness of S(T). 
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prototilc. The arclcngths 6 n f3 n r must go to zero in order for the prototilcs to be 
topological disks. 

Let d m ax and v m i n be the maximum diameter and minimum volume of prototiles 
in V, respectively. The tiles inside of S n (P) adjacent to the circular arc of length 
6(3 n r lie inside a strip of volume 

0/2[((3 n r) 2 - (J3 n r - d max f] + 2vd max 2 , 

where 27rdmax 2 ensures the ends of the strip are covered. The maximum number 
of tiles that can fit inside this strip is 

6 /2[2d max f3 v — dmax ] + 2,7rd max ^ ^d max f3 v ^ ^ 

Vmin V m in 

where K = (2n — 9/2)d 

max / v min does not depend on n. The minimum number 
of tiles required to cover the arc of length 9(3 n r is (8(3 n r)/(8 n f3 n r) = 8/9 n . For all 
n = 1, 2, ... we must have 

9 . 9d max f3 n r 
— < + K, or 

"n V m in 

9< d -^(6 n f3 n r) + K6 n . 

We have noted that both terms on the right must go to zero, forcing 8 to be zero 
also. This contradiction completes the proof. □ 

4.1. Local finiteness when the length expansion is Pisot. For the result 
we are going to prove, we will need some mild additional geometric hypotheses 
on prototiles and tilings. A straight edge I in the boundary of a tile T is a line 
segment that cannot be extended without leaving the boundary. Any such edge 
has a left and a right endpoint, relative to the boundary being traversed clockwise. 
The absence of circular arcs explains the terminology in the following definition. 

Definition 4.2. A tile T is called fractagonal if: 

(1) T has only finitely many straight edges, and 

(2) for any straight edge I ofT, the rest ofT lies in the interior of a half-plane 
generated by £ (we call this weak convexity,). 

The authors conjecture that the first property is automatically satisfied by the 
prototiles of a self-similar tiling substitution, provided they are weakly convex. 

We also need an assumption that for straight edges, collinearity implies comea- 
surability in the following sense: 

Definition 4.3. We say a tiling substitution S is a property-(C) substitution if 
whenever two straight edges l\ and £2 are collinear in a tiling t <E Xs, there is an 
r > for which \£i\, \£ 2 \ £ rQ(f3), where f3 = |A| is the length expansion of S. 

Theorem 4.4. Suppose that S is a property-(C) tiling substitution on a finite set 
V of fractagonal prototiles. If the length expansion of S is a Pisot number, then 
any tiling in r e Xs is locally finite. 

Notes. (1) We point out two common assumptions that are not necessary here. 
The substitution need not be primitive in that the structure matrix M may never 
satisfy M n > 0. Tilings in Xs need not be translationally finite, i.e. prototiles may 
appear in an infinite number of orientations. 
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(2) The authors are not aware of any primitive tiling substitutions that are not 
property- (C). 

(3) In his dissertation [3], Dirk Frettloh obtained a similar result under much 
stronger assumptions, including that the substitution is primitive, the expansion is 
an integer, and that the tiles are polygonal. 

Proof. A simple adjacency in t £ X$ is a two-tile patch a C r where the intersection 
is along straight edges. Kenyon ^ proves that for any tiling r made from a finite 
prototile set V, if t has an infinite number of inequivalent two-tile patches, then 
those patches are either simple adjacencies, or they occur along an entire circle of 
tile boundaries. By Lemma 14. II the latter cannot occur in a substitution tiling. It 
follows that in order to show r G X$ is locally finite, it suffices to show that there 
are only finitely many simple adjacencies a. 

Fix a simple adjacency a in r £ Xg- Two straight edges define a; the difference 
between two vertices, one selected from each edge, is called an offset of a. Thus 
each a defines four different offsets. It suffices to show that the total number of 
offsets is finite. 

Let the offset O of a be fixed. We want to see O as the difference of sums of 
edge lengths of tiles (see Figure EJl • As such, it is convenient to consider edges 
as being intervals on the positive real axis. We call a tile T placed below (resp. 
above) the real axis with a straight edge on the real axis a lower (resp. upper) tile. 
By applying a rigid motion, we can take edge representatives with their left (resp. 
right) endpoints at the origin. We call these edge sets £e and £ u . 

By weak convexity, the substitution S induces a one-dimensional substitution Se 
on Si in the sense of Subsection 12.11 It is clear that the length expansion (3 of S 
is the expansion of Se, so the one-dimensional substitution has a Pisot expansion 
constant. If T is the lower tile containing e £ £e, we simply rotate S(T) so that 
A • e is mapped to fJ ■ e and look at the lower tiles of (3 • e. Similarly we can create 
the one-dimensional substitution S u on £ u by considering the substitution on upper 
tiles. In particular, S u is just the reverse of Se. 

By property (C), the edge sets £e and £ u are partitioned into a finite number of 
subsets for which there are constants ri, ...,r% > with each rt/rj, i ^ j, irrational 
modulo (3, such that each edge length is in riQ(f3) for some i. Note that the 
application of Si or S u preserves these classes. Without loss of generality we may 
assume that the edges of the simple adjacency a have lengths in Q(f3) (i.e. r = 1). 

By the construction of r, there is some smallest n > for which there exists a 
tile T such that the simple adjacency a is a sub-patch of S n+1 (T). The intersection 
of the two tiles in a is contained in a line segment L that is composed of tile edges 
and is maximal in the sense that it cannot be extended without either entering the 
interior of some tile in S n+1 (T) or leaving S n+1 (T) entirely. We call L a fault line 
through the tile S n+1 (T). Since O depends only on a, we can use L to compute the 
value of O, even if the tiles surrounding L in S n+1 (T) do not appear in r (except 
of course those in a) . 

If n — 0, then S(T) can be moved to the positive cc-axis so that an endpoint of L 
lies at the origin. By the weak convexity assumption it follows that L can be seen 
as a union of edges of upper tiles and also as a union of edges of lower tiles. In this 
case L, along with the information about which level-one tile it is in and which tile 
edges define it, is called an initial segment. 
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If n > 0, then by the minimality of n, there must be two level-n tiles whose 
boundaries intersect a. From this we see that L must contain a line segment L' 
that is a union of edges of level-n tiles. Since the tiles are assumed to be fractagonal, 
weak convexity implies that if S n+1 (T) is moved so that L' in on the rc-axis, then 
L' is the union of edges of upper level-n tiles and also the union of edges of lower 
level-n tiles. Thus it is the image under S n of some initial segment L". 

Since there are only a finite number of level-one tiles up to equivalence, there 
are only a finite number of initial segments up to equivalence. Let £ be a set of 
representatives of these equivalence classes, restricted to those that are composed 
of edges with lengths in Q(/3), so that they are commensurate with the edges of 
a. Without loss of generality we choose these representatives to be on the positive 
real axis with an endpoint at the origin. For each L £ C, the endpoints of the 
lower tiles of L create a partition = ciq < a± < ... < a,k t — \L\, where kg is the 
number of lower tiles of L and \L\ is the length of L. Any such for any initial 
segment L will be called a prefix of L. Similarly, the endpoints of the upper tiles 
of L create a partition = bo < b\ < ... < — \L\ where k u is the number of 
upper tiles of L. A suffix of an initial segment L is \L\ — bi, where bi is in the 
partition of L. In the same manner, we may compute prefixes of Sg(e) for 
and suffixes of S u (e) for e £ S u , restricting again to edges with lengths in Q(j3). 
Since |5^(e)| = |/£u(e)| = P\e\, the prefixes will be partition elements of [0,/3|e|], 
and the suffixes will be j3\e\ minus partition elements of [0, /3|e|]. We collect all such 
prefixes into a prefix set Q px and such suffixes into a suffix set Q sx . It is important 
to note that both Q px and Q sx are finite. 

Since each P € V is fractagonal, there are only a finite number of straight edges. 
Thus there is a q G N so that all straight edge lengths in the class of the simple 
adjacency a are in |Z[/3]. This implies that the lengths of all initial segments, and 
the lengths of all prefixes and suffixes in Q px and Q sx are also in ^Z[/3]. Since 
offsets are the difference between a sum of upper edges and a sum of lower edges, 



We have already seen that O lies in a segment L' that is the image X n L" for 
some initial segment L". Let Lq be the representative in C of L". The tiles in 
a, which intersect L' = A™ • L", are equivalent to lower and upper tiles Tg and T u 
intersecting (3 n ■ L . In fact the edges of and T u are lower and upper edges that 
are edges in Sg(Lo) and S u (Lq) when we consider Lq in terms of its lower and upper 
tiles. The endpoints of edges in a used to compute O correspond to points a e Ti 
and b £ T u (see Figure EJ, so to compute O we need to compute \b — a\. Without 
loss of generality assume a < b. 



We may compute the left endpoint a in terms of powers of (3 applied to elements 
of Qpx] we do the same for the upper endpoint b using powers of applied to 
elements of Q sx . (See Figure HJ). Let p n £ Q px be the largest prefix of L for 



0c£z[fl. 








PT|L 



o 



Figure 5. The simple adjacency a is equivalent to T u U Tg. 
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which p n p n < a. If equality holds, we have computed a, otherwise we continue. 
We denote by e n -i the lower edge immediately following p n in Lq. This edge will, 
after application of SJ 1 , contain a (otherwise p n wasn't maximal). So there is a 
largest prefix of e n _i which we call p n ~i € Q P x for which n p n + /3 n f> n _i < a. If 
equality holds, we are done and the remainder of the prefixes will be the zero prefix. 
If not, we have that p n -\ is a prefix of <Sf(e„_i) for some edge e„_i £ £e, so the 
lower edge immediately following p n -\ in St(e n -\) will, after application of S n , 
contain a. We can continue in this fashion, choosing prefixes p n -2, ■■■,Po until we 
obtain that a = Y^k=o @ k Pk- 

Similarly we can compute the value of b. Let s n € Q sx be the largest suffix of 
Lq for which b < f3 n (\Lo\ — s n ). If equality holds we are done and all subsequent 
suffixes are the zero suffix. If not we select the edge e n -\ immediately preceding 
|i | — s„ and find its suffix s n -i so that b < f3 n (\L \ — s n ) — /3™~ 1 s„_i. We continue 
choosing suffixes s„_2, so until we obtain that b = f3 n \Lo\ — X)fc=o P ks k- 



Ji-1 



a 



b 



„n. 

PlL ( 



Figure 6. The computation of a and b as seen in [0, /3™|£o|]- 



Thus the offset distance O is given by 

n 

(10) O = (3 n \L \+Y / P k (sk-p k ). 

fc=0 

Our interest, then, is in all possible combinations of suffixes and prefixes. We collect 
these into a difference set V = {— s — p such that s £ Q SX ,P S Qpx} which we must 
show is finite. 

If /3 is an integer, then the result follows immediately since O is bounded by the 
length of the longest straight edge and is an element of |Z. 

If {3 is not an integer, we need to use the ideas and notation of Section|21 Recall 
that the map (j> : C d — > C is a bijection between Q d and Q(/3), and that multiplica- 
tion by the matrix A in Q d is equivalent to multiplication by (3 in Q(/3). We can 
find preimages Wk € Q d such that <p(wk) = —Sk —Pk for k = 0, n. Since elements 
of T> are in -Z[0], we have that each Wk must be in Equation (|f Oj l becomes 

n / n \ 

(11) O = /?"|Lo| + X]^(^) = / 3n l £ o| + ■ 

fc=0 \k=0 J 

We write the preimages of T> restricted to Q d in terms of the eigenvectors 
v\, ...,Vd- Since <fi projects onto vi sending all other eigenvectors to 0, and since O 
is necessarily less than the length of the longest edge, there is an mi > that is an 
upper bound on the modulus of the component of _1 (O) fl Q d in the v\ direction. 
For any j — 2, d, denote by rrij the maximum modulus of the component in the 
Vj direction of any preimage of the difference set V restricted to Q d . The modulus 
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of the component of X)fc=o ■A^Wk m the Vj direction is less than or equal to 



n oo 



fc=o fe=o 1 

since |/3j| < 1 by assumption when j = 2, ...,d. So we see that the vectors which 
project onto O must lie in a bounded region. There are only a finite number of 
elements of in a bounded region, so there are only a finite number of offsets 
0, and hence only a finite number of simple adjacencies A. The tiling x is locally 
finite. □ 

Comment. A complex algebraic integer A is called a complex Pisot number if 
|A| > 1 and all algebraic conjugates A' ^ A satisfy |A'| < 1 (see jH], |E|)- If A is 
complex Pisot, (3 — |A| may not be real Pisot (see for an example). We do 
not know whether Theorem 14.41 always holds in the case of a complex Pisot linear 
expansion A. 

4.2. The Figure 111 example is not locally finite. Let a denote the large square 
prototile from Figured Since the substitution of any tile includes a tile of type a, 
it follows that a block of type S n (a) must appear in S' n+1 (e) for any tile type and 
for any n = 1, 2, ... and so every tiling x G X must contain blocks of type S n (a) for 
all n = 1,2, .... 

The circled point in Figure|5]appears inside the block S(a). This point represents 
a vertex of the tiles above it and it intersects the upper boundary of the tile below 
it at a point [3—1 units from the left endpoint of that boundary. Each time 
we apply the substitution, the circled point determines the offset between the tile 
below it and the two tiles above it. This offset is completely determined by F. If 
we can show that the orbit of (3 — 1 is infinite then this implies that there are an 
infinite number of offsets and hence the tiling is not locally finite. It is convenient 
to compute this orbit for / rather than F. 

Example gives the algebraic generalized /3-transformation associated with the 
substitution on upper edges of tiles. Noting that the eigenvalues are real, we com- 
pute the corresponding action / in We C R 2 . The companion matrix of (3 is 



A 



1 

1 3 



We choose as the generator of V + the vector v = (—(32, 1) and as generator of Vq 
we take V2 = {—f3, 1). The vectors mapping onto 2/1,2/2,2/3,2/4 and 2/5 are given by 
m = (0,0), w 2 = (0,0), w 3 = (1,0), w 4 = (2,0), and w 5 = (3,1). The sets W 3 
have boundaries given by Vq + 4>^ 1 (xj-i) and Vq + (f>^ 1 (xj), where the inverse is 
taken in Q 2 . It is not difficult to check that the appropriate translation vectors are 
(0,0), (1,0), (2/3, 1/3), (1/3, 2/3), (0,1) and (1,1). (In this example one can also 
find these by looking at A~ 1 Wj.) The regions Wj are as shown in Figured along 
with the main eigenspace V + . 

We now establish an /-orbit which enters the escape region and thus tends to 
infinity. By our theorem, this orbit corresponds to an infinite f-orbit. We begin 
by computing the escape region. Looking in the proof of Theorem 13.11 , we see 
that we must choose T2 so that 

r2> LTT^T fora11 J = 1:2,3,4. 
IP2I - 1 
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FIGURE 7. The regions Wj in the domain W E of /. 



We recall that w\ = (0,0), «J 2 = (1,0), W3 = (2,0), and w 4 = (3,1) and compute 
their components in the Vq direction: 



W1.2 — 0, W 2 ,2 = 



■1 



W 3 ,2 



I0 4 o 



-7 + Vl3 



T3 "'" V13 ' 2V13 
The largest of these in modulus is \ws^\ ~ -5547. We see that 



r 2 > 



Iftl- 



1.8321, 



so we simply choose r 2 = 2 to define lZ esc . 

We will show that w = (p~ 1 (f3 — 1) = (—1,1) enters the escape region by com- 
puting the first several elements of its orbit. (This was made easier by using a 
cobwebbing program in Matlab to determine in which Wj the point was at each 
step.) Now (-1, 1) e Wi, so /(-l, 1) = A(-l, 1) - Wi = (3,0). One can compute 
that (3, 0) € W4, so /(3, 0) = A(3, 0) — W4 = (—3, 2). Continuing in this fashion the 
next several iterates are (5,-1), (—6, 3), (9, —3), (—11, 6), (15, —6). The orbit enters 
the escape region at (—6,3) and the iterates begin to grow quite quickly. It is in- 
structive and surprising to draw the region We, which is infinitely long and quite 
thin, and plot these points which must remain inside but which bounce further and 
further away from the origin. 
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